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Abstract. Integrability, algebraic structures and orthogonal basis of the Calogero model 
are studied by the quantum Lax and Dunkl operator formulations. The commutator al- 
gebra among operators including conserved operators and creation-annihilation operators 
has the structure of the W-algebra. Through an algebraic construction of the simultane- 
ous eigenf unctions of all the commuting conserved operators, we show that the Hi- Jack 
(hidden- Jack) polynomials, which are an multi-variable generalization of the Hermite 
polynomials, form the orthogonal basis. 



1 Introduction 



In memory of the pioneering works in 70 's |0, ^ul| , |Su2|| , a class of one- dimensional 



quantum many-body systems with inverse-square long-range interactions are generally 
called the Calogero-Moser-Sutherland models. The celebrated Hamiltonians are 

1 ^ 1 ^ - ha 
Calogero-Moser: i^cM = 77 E^i + o E 7 yi^ (1) 

Calogero: fl, = \ ^(p? + + \ ^ (2) 

^ i=l ^ i,fc=l ^-^3 -^k) 

jj^k 

Sutherland: = J ^^pj + ^ —4? (3) 

)+k 

where the constants iV, a and lo are the particle number, the coupling parameter and 
the strength of the external harmonic well, respectively. The momentum operator is 
given by a differential operator, = —\h-^. The Calogero and Sutherland models are a 
harmonic confinement and a periodic version of the Calogero-Moser model, respectively. 
Thus these two models have discrete energy spectra, whereas the other has continuous 
one. From now on, we set the Planck constant at unity, h = 1. 



Based on the talk given by MW at the Workshop on the Calogero-Moser-Sutherland models at CRM, 
Montreal, March, 1997. To appear in the proceedings of the workshop. 
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The Lax formulation for the classical Calogero-Moser model was discovered by Moser ||M 
Let us introduce two N x N Hermitian matrices: 

Lij = PiSij + ia(l - 6ij) — - — , 

N ^ ^ 

We call them Lax pair. The classical Calogero-Moser Hamiltonian is given by eq. ([l|) 

dXj_ 

dt 



with pj = ^ and h = 0. The time evolution of the L-matrix is expressed as the Lax 
equation, 

^ = {L,H^^^}^ = [L,iM], (4) 

where the Poisson bracket is defined by {/, g}p =^ Thanks to the trace 

identity for c-number-valued matrices, TiAB = TtBA, we can easily see that the trace of 
the power of the L-matrix, '^TiL"', gives the conserved quantities: = Tr[L", iM] 



0. Using the classical r-matrix [^] or the generalized Lax equations for higher conserved 



quantities ||BK|| , we can show that the conserved quantities are Poisson-commutative, 
{/^',/^}p = 0, which proves the integrability of the classical Calogero-Moser model in 
Liouville's sense. Natural quantization of the Lax equation for the classical case (^) 
by correspondence principle, <|fc}p — > — i[J!k, 4], gives an equality for the quantum 
Calogero-Moser Hamiltonian (^. However, the trace trick is not available to construct the 
commuting conserved operators for the quantum case because of the non-commutativity 
of the canonical conjugate variables. The initial motivation of our study was to find 
out a way to construct the conserved operators for the quantum models using the Lax 
formulation. The key of our idea is the sum-to-zero condition of the M-matrix: 

N N 

Y,M,k = 0, Y.Mkj = 0, forfc = l,2,---,iV. 

j=i i=i 

This property tells us that the (commuting) conserved operators can be obtained by 
summing up all the matrix elements of the powers of the L-matrix instead of taking 



traces UHW, UWH 



N 

= Y: = TsL" ^ [HcM,!^''] = Ts[L",M] = 0, 

j,k=i 

which proves the quantum integrability, or the existence of sufficiently many conserved 
operators, of the Calogero-Moser model. Encouraged by the result, we further investigated 
the integrable structure of the Calogero model (0) through the quantum Lax formulation. 

2 Integrability and Algebraic Structure 

Let us start from the Lax equation of the Calogero Hamiltonian (^), 
-i^ = [Hc,L^] = [L^,M]±u;L^, 
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where the new matrices L'^ are defined by L^'^= L±. Q, Qjk'= ixjdjk- Using the sum-to- 
zero trick, we can get the conserved operators of the Calogero model as follows: 

L = Tj,{L+L-r, [HcJn] = Ts[(L+L-)",M] = 0, 



h = 2Hc - Nu;{Na + (1 - a)), /„ = J^Pf + 



N 



Mutual commutativity of the above conserved operators is verified rather easily by the 
Dunkl operator formulation [0, Introducing the coordinate exchange operator, 

iKikf){- ■ ■ ,xi, - ■ ■ ,Xk, - ■ ■) = f{- ■ ■ ,Xk, - ■ ■ ,xi, - ■ 
we define the creation-annihilation like operators as 

N ^ 

cj = pi + iaY^ Kik + iujxi, 



k=l 

N 



ci = pi + iaV] Kik - iujxi. (5) 

Commutation relations among the creation-annihilation operators are 

N 

[ci, clJ = 2uj6im{l + aYl ^ik) - 2u;a(l - 6im)Kim, 

k^l 
k=l 

which prove that the Hermitian operators. In *== Z]jLi(c]cj)", are commuting operators. 
We denote the restriction of the operand to the space of symmetric functions by ^ . 

Under the restriction, the conserved operators J„ and commuting Hermitian operators 
l„ are considered to be the same, /„ = l„ • Thus we have proved the quantum 

Sym Sym 

integrability of the Calogero model. 

We can recursively construct generalized Lax equations for a family of operators O^, 
m,p = 1, 2, ■ ■ -, which reveal the W-algebraic structure of the Calogero model | UW1 |. The 
operators are defined by the sum of all the matrix elements of the Weyl ordered product 
of p L+s and m L~s: 



Jw 



[(L+)^(L-r]w = 7-^ E {LmL-r- 

' ' all possible 
order 

The generalized Lax equations for the operators are 

[Oj;,L±] = [L±,M^]+mc.(l±l)[(L+nL-)™-V 



-pu{\^\)\[L^r-\L-r\^, (6) 
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and the matrices satisfy the sum-to-zero condition: 

N N 

i=i i=i 

The Hamiltonian Hq belongs to the operator family, 2Hc = 0|. The operators O" 
are conserved operators, though they do not commute each other. The family has two 
interesting subsets of commuting non-Hermitian operators, 

i?r=Oo", Bj^Ol n = l,2,---, (7) 

which we call power-sum creation-annihilation operators. They will play an important 
role in the algebraic construction of the energy eigenf unctions in the next section. 
Let us introduce operators Wjf^: 

Wi^^'='^0:-:zl s>\n\ + l, 

where the indices n and s are integer or half odd integer, and respectively correspond to 
the Laurent mode and the conformal spin. The commutator among the operators above 
is 

[w^i^^ w^i!^] = Ht - 1) - mis - i))w!iZ''^ + vi:'^iwt\ 

where 'P^^i';^{Wi^^) is a polynomial of Wi^\ u<s + t — 3,l<n + m. The polynomial 
is generated while the products of L^-matrices are rearranged into the Weyl ordered 
products by replacements of L"*" and L~ , 

L~] = 2uj{{a - 1)1 - aT), Ijk = 5jk, Tjk = 1. 

In terms of the W-operators, conserved operators and power-sum creation-annihilation 
operators are respectively expressed as oc Wq~^^, oc W_l, and B^ oc Wn^ 

For the classical case, the W-algebraic structure of the Calogero model was discovered 
by the collective field theory ||AJ 1|| and the classical r-matrix method . The quantum 



collective field theory also possesses the W-algebraic structure [|AJ2| , |A J3|| , though its 



relationship with the quantum Calogero model is not directly confirmed. An SU{iy) 



generalization of our approach is presented in [LJW2] 



3 Perelomov Basis 



The eigenvalue problem of the Calogero model was first solved by Calogero . Later, 
inspired by the simple form of its energy spectrum, Perelomov tried an algebraic construc- 
tion of the energy eigenfunctions |[Pe|| . In what follows, we shall complete the Perelomov's 
approach ||U Wl|| . The generalized Lax equations (|^) for the power-sum creation operators 



yield the following commutators: 
[Hc,Bi]=nuBl [BlBl]={ 



n, m 



1,2, 



N. 
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To construct all the eigenf unctions for the A^-body Calogero model, we need creation 
operators. By straightforward calculations of the commutators (§), Perelomov presented 
three creation operators, Bj^, n = 2,3, 4. The quantum Lax formulation provides an easy 
way to make sufficient number of such operators. By successive operations of the power- 
sum creation operators on the ground state wave function, we can get all the excited 
states. The Calogero Hamiltonian is cast into the following form. 



1 1 1 ^ 

= 2 Te L+L- + -Nu{Na + (1 - a)) = - ^ h]h, + E, 



where the operators and hj are defined as 

N ^ I 

h] -i? Ltj = Pj + iujXj -iaJ2 - — 



k=l k=l '^i 

N ^ 1 



fc=i k=i i ^ 

k^j 

Thus the differential equations, hj\Q) = 0, j = 1,2, ■ ■ ■ , N, are the sufficient conditions 
for the ground state |0). In the coordinate representation, the solution is expressed as the 
(real) Laughlin wave function: 

1 ^ 

(x|0)= n |x,-Xfcrexp(--a;E^|)- (9) 

^<j<k<N ^ j=l 

As is similar to the free boson case, the excited states are labeled by the Young 
diagram. 



A = {Ai > As > ■ ■ ■ > Ajv > 0} G Y/ 



N, 

where A^, k = 1,2, ■ ■ ■ , N, are nonnegative integers. Conventionally, we omit zeroes and 
use superscript for a sequence of the same numbers, e.g., {4, 1^} = {4, 1, 1, 0, ■ ■ ■ , 0}. The 
excited state labeled by the Young diagram A is given by 

N 

|A) = 11(^1)'^"''=^^ |0), A;v+i=0, 



k=l 



Hc = i\X\uj + E,)\X)''^'EiX)\X), (10) 

where |A| denotes the weight of the Young diagram, \X\'^ J2k=i ^k- The above energy 
spectrum has the same form as that of non-interacting bosons confined in an external 
harmonic well up to the ground state energy. In other words, the inverse- square interac- 
tions just shift the ground state energy. The multiplicity of the n-th energy level, nu + Eg, 
is equal to the number of the Young diagrams of the weight n, #{A||A| = n}. This is also 
the same as that of the non-interacting case. Thus we have algebraically constructed a 
basis of the eigenfunctions of the Calogero Hamiltonian. 
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4 Diagonalization of I2 



The algebraic construction a la Perelomov generates a basis of the eigenf unctions 
of the Calogero Hamiltonian. Unfortunately, the basis is not orthogonal. To make an 
orthogonal basis from a basis, we usually try the Gram-Schmidt method. Here we take 
another way. As we have confirmed before, the Calogero model has a set of commuting 
conserved operators, which means existence of simultaneous eigenfunctions for them. The 
simultaneous eigenfunctions should be the orthogonal basis because they must be non- 
degenerate eigenfunctions of Hermitian operators. As the first step of our approach, we get 
some simultaneous eigenfunctions of the Hamiltonian and the second conserved operator 



I21 and observe their properties |PJW3 | 



Since the Hamiltonian and I2 commute, the matrix representation of I2 on the Perelo- 
mov basis has a block-diagonalized form and each block consists of the wave functions of 
a weight (energy eigenvalue). By a straightforward calculation of commutators between 
I2 and we calculated the first seven blocks, whose weights are from zero to six, of the 
matrix representation and their eigenvalues. The eigenvalues imply the general form of 
the eigenvalue of J2: 

N 

E2{\) = Au^ Y.({\kf + a{N + l- 2k)\k). 

k=l 

Though a combination of E{X) and -E'2(A) removes most of the degeneracies, there still 
remain degeneracies for the states whose weights larger than or equal to six. For example, 
the following two pairs of the Young diagrams with the weight six give such degeneracies: 

{4,1^}, {3^} ^ E2 = Au;\l8 + 6a{N -2)), 
{3,1^}, {2-'^} ^ ^2 = 4cu=^(12 + 6a(A^-3)). 

It is interesting that the pairs have a common property. We can not compare two Young 

D 

diagrams of each pair by the dominance order. The dominance order <, sometimes called 
the natural partial order, is defined as follows: 

D ' ' 

fj,<\ -v^ \fj.\ = \X\ and ^ /ifc < ^ Ayfc for all /. 

k=l k=l 

We can readily confirm that the Young diagrams of each pair are incomparable in the 
dominance order, 

{4,in|{3nand {3^ § {4, 1^, 
{3,1^} I {2^} and {2^} § {3,1=^}. 

The specific observation above is, in fact, a general fact. We cannot define the dominance 
order between any pair of distinct Young diagrams A and /i of a weight that share the 
common eigenvalue E2 P^, |U W6||, i.e.. 



D D 

|A| = and -E2(A) = E2{fi) ^ X fi and /i ^ A. 
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We calculated the eigenvectors of the blocks with weights up to three in the matrix 
representation of J2. The eigenvectors correspond to seven simultaneous eigenf unctions of 
Hq and l2- Since the eigenvalues E and E2 for the seven functions have no degeneracy, 
they belong to the orthogonal basis and also to the simultaneous eigenfunctions of all 
the commuting conserved operators of the Calogero model. The eigenfunction of the 
Calogero model is factorized into the ground state wave function @ and a symmetric 
polynomial. Symmetric polynomial parts of the seven simultaneous eigenfunctions, which 
we denote by [A], are 



2, , a N{N - 1) 



[0] = 1, [l]=mi, [l^]=mi. + 2^ 2 , 
[2] = (1 + a)m2 + 2ami2 - —N{Na + 1), 
[A = rni3 + —a^ ^mi, 

[2, 1] = (2a + l)m2,i + Garriis - ^{1 - a){N - l){Na + l)mi, 

[3] = (a^ + 3a + 2)m3 + 3a(a + l)m2,i + Garriis 

- —(a^N^ + 3aN + 2)mi, 
2uj 

where mx is the monomial symmetric polynomial defined by 



ctSSat, distinct 
permutations 



Note that the sum over Sn is performed so that any monomial in the summand appears 
only once. In the above expressions, we notice that all the seven symmetric polynomials 
share a common property, triangularity. Namely, the seven polynomials [A] are expanded 
by the monomial symmetric polynomials whose Young diagram /i is smaller than or 

d 

equal to the Young diagram A in the weak dominance order <, i.e., 
/i<A ^22l^k<22'^k for all /. 

k=l k=l 

The observation means that we can uniquely identify the simultaneous eigenfunctions of 
the first two conserved operators of the Calogero model just by the first two eigenvalues 
and triangularity up to normalization. We shall confirm the existence of such functions 
by algebraically constructing them. 



5 Hi-Jack Polynomials 

Since our interest now concentrates on the symmetric polynomial parts of the simul- 
taneous eigenfunctions, we modify some operators to make them suitable for the aim. A 
gauge transformation of the creation-annihilation-like operators (§) yields the following 
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Dunkl operators: 



t def 

a = 



i ^ 1 

— {pi + iaY^ {Kik - 1) + 2iuxi) 



2^^ xi - Xk 



k^l 

1 ^1 

ai = (x|0)iQ— — - = i{pi + ia 51 i^ik - 1)), 

{x\0} xi - Xk 

k^l 

di = a]ai. (11) 

The gauge transformation above removes the action on the ground state wave function 
from the operators. Note that the definition of Hermiticity of such gauge-transformed 
operators is modified and different from the ordinary one. Using the ^/-operators, we 
define the normahzed conserved operators: 

J^l ^^"^ ^2CJ^ ^ ' ' (x|0) Sym' 

(^|0)^C^=^/l+i?g. (12) 

We note that the Dunkl operators (pJ]) reduce to those for the Sutherland model (|^) in 
the limit, uj ^ oo: 

a\ zi, 

N I dci d 

^1 - Zk dzi 

k^i 

di^Di = ziVi. (13) 

We change the variables by 

exp2ixj = Zj, j = 1,2, ■ ■■,A^, 

and denote the ground state wave function and the ground state energy of the Sutherland 
model by 

N ^ 

4 = n \z,-zk\'^xiz-^^'^-'\ 

^<j<k<N j=l 

eg = ^a^N{N~l){N+l). 

Then the Sutherland Hamiltonian (^ is gauge-transformed to and related with the D- 
operator by 
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j,k=i "J 



Zj Zj^ 



N 



2E(A 



1=1 



ISym 



Commutation relations among the Dunkl operators (0) and the action of ai on 1 are 



N 



k=l 



[di, dm] = a{drn - di)Kim, a; ■ 1 = 0. (14) 

We should remark that the above relations do not explicitly depend on the parameter 
a;, which implies the Dunkl operators for the Sutherland model (|l^) also satisfy the 
above relations. Hence the Calogero and Sutherland models share the same algebraic 
structure [[UW4| , |UW5| , [Ul|| . To put it another way, the theory of the Calogero model 
is a one-parameter deformation of that of the Sutherland model. Thus the simultaneous 
eigenfunction of the Calogero model is expected to be a one-parameter deformation of 
that of the Sutherland model, which is known to be the Jack polynomial In the 
following, we call the simultaneous eigenfunction of the Calogero model Hi- Jack (hidden- 
Jack) polynomial. 

Using the normalized conserved operators (0), we define the Hi- Jack polynomials 
jx{x] tu, 1/a) in a similar fashion to a definition of the Jack polynomials: 

N 

Iiix{x]UJ,l/a) = J2^kjxix;uj,l/a) = Ei{X)jxix;uj,l/a), 



(15) 



k=l 
N 



l2jx{x;uj,l/a) = Y.{Xl + a{N + l-2k)Xk)jx{x;uj,l/a) 



k=l 



def 



E2{X)jx{x;uj, 1/a), 



(16) 
(17) 



d 



1. 



wxxia,uj) 

We can prove the existence of the Hi- Jack polynomials by explicit construction. Fol- 
lowing the results by Lapointe and Vinet on the Jack polynomials [CV]], we introduce the 
raising operators for the Hi- Jack polynomials, 

K= E a%^j, fork = 1,2,- ■■,N-l, 

JC{1,2,--,N} 
\J\=k 



The operators, a\ and dij, stand for 



a : 



a 



J' 



d 



i,j 



n 

{dj-^ + a) {dj2 + 2a) ■ ■ ■ {dj^ + ka) 
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where J is a subset of a set {1, 2, ■ ■ ■ , A^} whose number of elements \J\ is equal to k, 
J C {1, 2, ■ ■ • , N}, \ J\ = k. From eq. (|14|), we can verify an identity, 



{di + ma) {dj + {m + 1] 



Sym 



{dj + ma) {di + {m + 1) 



Sym ' 



(19) 



where m is some integer. The symbol where J is some set of integers means that the 

Sym 

operands are restricted to the space which is symmetric with respect to the exchanges 



of any indices in the set J. This identity (JT^) guarantees that the operator dij does 
not depend on the order of the elements of a set J when it acts on symmetric functions 
and hence operation of the raising operators on symmetric functions yields symmetric 
functions. The function generated by the following Rodrigues formula, 



Jx{x;uj,l/a) = C^\b+)'-{b+,,)'--^-'- 
with the normalization constant Cx given by 



Ai— A2 



Af-1 



ca — n ^^(-^1' -^2 



Afc+i; a), 



k=l 



where 



A2, ■ ■ ■ , Afc+i; a) — (a)Aj.-Afe+i(2a + Afc_i — Afc)Afe-Afc+i 
■■■(A;a + Ai-Afc)A,-Afe+i, 

satisfies the definition of the Hi- Jack polynomial j\{x;uj,l/a). The symbol (/?)„ in the 
above expression is the Pochhammer symbol, that is, {l3)n = P{P + I) ■ ■ ■ {P + n — 1), 

{P)o = l. 

The first seven Hi- Jack polynomials are, for instance, given as follows: 

jo{x;uj, l/a) = Jo{x; 1/a) = mo{x) = 1, 
Ji{x;uj, l/a) = Ji{x; l/a) = mi{x), 

ji2{x;uj,l/a) = Ji2{x; l/a) + ^-^^ — —Jo{^'^ V^) 

, , a N{N-1) . . 
= mi2{x) + — ^ mo(a;), 

{a + l)32{x\ l/a) = {a + l)J2{x; l/a) - ^N{Na + l)Jo{x; l/a) 

= (a + l)m2{x) + 2ami2{x) N{Na + l)mo(a3), 

1 {N -1){N -2) , 
ji3{x;u, l/a) = Ji3(x; l/a) + —a Ji{x; l/a) 

. . 1 {N-l){N-2) , , 
= m^z[x) H a m\[x), 

{2a + l)32,\{x\i^,l/a) 

= {2a + l)J2,i{x; l/a) - —(1 - a){N - l){Na + 1) Ji(a;; l/a) 

2uj 

= {2a + l)m2,i(a;) + 6ami3{x) 
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- ^(1 -a){N- l)(iVa + l)mi(a^), 

{a^ + 3a + 2)j3{x; lu, 1/a) 

= (a^ + 3a + 2)J3{x; 1/a) - —{a^N^ + 3aN + 2)Ji{x; 1/a) 

2uj 

= (a^ + 3a + 2)7713(2;) + 3a(a + I)m2,i(a3) + Ga'^niis^x) 
3 

- — (a^iV^ + 3aA^ + 2)mi{x), 
2uj 

where the symbol J\{x] 1/a) denotes the Jack polynomial. The explicit forms also show 
the fact that the Hi- Jack polynomial reduces to the Jack polynomial in the limit, u ^ 00, 

jxix;iu = 00,1/a) = Jxix;l/a). (20) 

Besides the above relation, we have some other relations between the Hi- Jack and Jack 
polynomials [|U W5| , |D1|, |U2|| . While the Hi- Jack polynomial is a one-parameter deforma- 
tion of the Jack polynomial, we can get the Hi- Jack polynomial from the Jack polynomial 
by the following formula, 

Jx{ai,(yl, ■ ■ ■ ,a]v; 1/a) ■ 1 = Jx{x;uj, 1/a), 

which gives another relation between the Jack polynomials and the Hi- Jack polynomials. 
In the above expansion, we have an observation, and it is generally true, that increasing 
the order of 1/2uj by one causes decreasing of the weight of the symmetrized monomial 
by two. The fact yields an stronger form of the triangularity: 

jx{x;uj,l/a) = 2^ i^—j wx^^{a)mf,{x), 



d 

IJ,<X and 
|At| = |A| (mod2) 



wxxia) = 1. (21) 



Combining eqs. ( PDj ) and (PT|), we have the following expansion form of the Hi- Jack 
polynomial with respect to the Jack polynomial: 

Jxix;uj, 1/a) 

= Jx{x- u, 1/a) + (^)^'"-'^'^/W(«)^m(^; !/«)• 

d 

li<\ and |/i|<|A| 
and |/i| = |A| (mod2) 

Relationship between the Hi- Jack polynomials and the Perelomov basis (|iy) is given as 
follows. The power-sum creation operator is cast into the power-sum of aj-operators. 



N 

Bi = i2iu;)'^Yi^l)± =(2io;)'=p.(at) 

Sym 



Sym' 



1=1 

and the Perelomov basis is expressed by the power-sum of the aj-operators as 
^^'^^ - (2i.;)l^l n(Pfc(at))^-A^-^i . 1 = (2i^)l^lp,(at) . 1. 



(^|0) 



fc=i 
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Thus the transition matrix between the power-sums and Jack polynomials M(J, p), 

J\{x,l/a)= M{J,p)x^,Pf,{x), 

ImI=|a| 

gives a relation between the Hi- Jack polynomials and the Perelomov basis, 
Ux;uj,l/a) = {2iu;)-^'\ J] M{J,p)J^. 

H = |A| 

We have introduced the Hi- Jack polynomials as the simultaneous eigenfunctions for the 
first two commuting conserved operators with the triangularity. As we shall see shortly, 
they are non-degenerate simultaneous eigenfunctions for all the commuting conserved 
operators of the Calogero model. From a calculation of the action of di operator on a 
symmetrized monomial of a|,'s, mx{ai, ■ ■ ■ , ajv), we can prove the following expression: 

Injx{x;uj,l/a) = J2 ^A,M(«'l/2t^)"^/^(*)- (22) 

D 

or |/^|<|A| 

This means that operation of the conserved operators on the Hi- Jack polynomials keeps 
their triangularity. Since the n-th conserved operator commutes with the first and second 
conserved operators, [/i,/n] = [h,In] = 0, we can easily verify, 

Iilnjx{x;uj,l/a) = Ei{\)Injx{x;uj,l/a), (23) 
hinjxix; uj, 1/a) = E2{X)Injx{x; uj, 1/a). (24) 

Equations (0), (p^) and (p2|) for Injx are respectively the same as eqs. (pTSf), (|16D and 
(|l^) for the Hi- Jack polynomial j\, which means Injx satisfies the definition of the Hi- Jack 
polynomial except for normalization. Our definition of the Hi- Jack polynomial uniquely 
specifies the Hi- Jack polynomial. So we conclude that Injx must coincide with jx up to 
normalization. Thus we confirm that the Hi- Jack polynomials jx simultaneously diago- 
nalize all the commuting conserved operators /„, n = 1, - ■ ■ ,N. The eigenvalues of the 
conserved operators, 

Injx{x;uj,l/a) = En{\)jx{x;uj,l/a), 

are generally polynomials of the coupling parameter a: 

^n(a)=ef(A) + e«(A)a + ---. 

It is easy to get the constant (a-independent) term e^\\) because the term corresponds 
to the n-th eigenvalue for non-interacting bosons confined in an external harmonic well: 

N 

er(A) = E(^^)'^- 

k=l 
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It is clear that there is no degeneracy in the constant terms of the eigenvalues {e^\X)\n = 
1, ■ ■ ■ , A^}. Since the conserved operators J„ are Hermitian operators concerning the inner 
product, 



{jxJ/i) = n ^^k\{x\0)\'^jxjf^ oc 5x,^, 



the Hi-Jack polynomials are the orthogonal symmetric polynomials with respect to the 
above inner product. From the explicit form of the weight function, 



N 

2 



l(^|0)P= n \xj-Xk\'''exp{~ujY.x', 

^<j<k<N 1=1 

we conclude that the Hi- Jack polynomial is a multivariable generalization of the Hermite 
polynomial ||La| . 
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